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ABSTRACT 

An abstract approximation framework for the identification of nonlinear distributed parameter 
systems is developed. Inverse problems for nonlinear systems governed by strongly maximal mono- 
tone operators (satisfying a mild continuous dependence condition with respect to the unknown 
parameters to be identified) are treated. Convergence of Galerkin approximations and the corre- 
sponding solutions of finite dimensional approximating identification problems to a solution of the 
original infinite dimensional identification problem is demonstrated using the theory of nonlinear 
evolution systems and a nonlinear analog of the Trotter-Kato approximation result for semigroups 
of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing 
linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic 
operators and the corresponding nonlinear parabolic partial differential equations to which they 
lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is 
discussed. 
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1. Introduction 

In this paper we develop a general abstract approximation framework 
for the identification of nonlinear distributed parameter evolution systems. 
Our intent is to define relatively straightforward and easily verified 
criteria that are applicable to broad classes of nonlinear systems; these 
criteria will guarantee the convergence of solutions to a sequence of finite 
dimensional Galerkin approximation based parameter estimation problems to 
a solution of the original, underlying, infinite dimensional identification 
problem. The results that we present below generalize and extend the 
theory recently developed by Banks and Ito in [2] and [3] for regularly 
dissipative or abstract parabolic, linear systems. It is, to the best of our 
knowledge, the first such general approximation theory for inverse 
problems involving nonlinear distributed systems. 

The sufficient conditions set down in our framework include a 
relatively mild continuity assumption with respect to the unknown 
parameters to be identified, an equi-boundedness and an equi-strong 
monotonicity assumption on the nonlinear operator describing the system 
dynamics. In addition our theory requires a standard approximation 
assumption on the Galerkin subspaces used to effect the finite dimensional, 
or finite element, approximations. We demonstrate that solutions to the 
finite dimensional identification problems approximate a solution to the 
infinite dimensional identification problem via a convergence result for 
solutions to the forward problems. This result is obtained using the theory 
of nonlinear evolution systems and a nonlinear analog of the well-known 
Trotter- Kato approximation result for linear semigroups. 
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In the present paper, we are concerned only with theory; 
implementation questions and conclusions drawn from our numerical or 
computational studies will be reported on elsewhere. Also, while we have 
tried to make our framework as versitile as possible, the treatment below 
does have limitations. For example, our theory can handle quasi-autonomous 
systems but it is not applicable in the fully nonautonomous case. The 
development of a general theory which can handle nonlinear systems 
involving time dependent operators requires additional effort and is 
currently the focus of our ongoing investigations. The particular 
difficulties inherent in the time dependent case will be described in 
greater detail in our discussions below. 

We provide a brief outline of the remainder of the paper. In Section 
2 we state a fundamental existence and uniqueness result for infinite 
dimensional nonlinear systems and prove a general approximation result 
which is especially well suited for application in the context of the inverse 
problems which are the central focus of our study. In Section 3 we 
define a class of nonlinear distributed systems and the associated 
parameter identification problems. We define the Galerkin approximations 
and prove the general convergence result. Section 4 contains some examples. 
We show that our nonlinear theory subsumes the linear theory presented in 
[2] and [3] as a special case; we also consider the application of our 
framework to a class of nonlinear elliptic operators and the corresponding 
nonlinear parabolic partial differential equations to which they lead. In 
particular, we look at the application of our results to a well known 
quasilinear model for heat conduction or mass transfer. In Section 5 we 
summarize our findings and provide some concluding remarks. 
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2. An Approximation Result for Nonlinear Evolution Systems 

Let X Q be a Banach space with norm || 0 . We consider the nonlinear, 
quasiautonomous initial value problem in X Q given by 

(2.1) x 0 (t) + A 0 x 0 (t) 9 f 0 (t), 0 < t tT, 

(2.2) x Q (0) = Xq 

x 

where x“ € X 0 , f Q : [0,T] -» X Q and the nonlinear operator A Q : X Q -* 2 0 is in 
general multivalued, not everywhere defined, and not continuous. The 
existence of solutions to the initial value problem (2.1), (2.2) and the 
subsequent approximation result to follow, are both consequences of Theorem 
2.1 to be given below. 

We shall require the following definitions. Let X be a Banach space with 
norm | | x . For A: X -» 2 X , a nonlinear, multivalued operator, the domain and 
range of A are defined by Dom(A) = {x € X: Ax / 0} and R(A) = u Ax 

X6dom(A) 

respectively. We say that the operator A is accretive if for every X > 0, 
Xj,x 2 e Dom(A) and y x € Ax r y 2 € Ax 2 we have 


i x i - x 2'x * i x i ~ x 2 + x (yi “ y 2 )i x - 

We say that A is m-accretive if A is accretive and K(I + XA) = X for some 
X > 0. We note that if A is m-accretive then R(I + XA) = X for every X > 0 
and for each X > 0 the resolvent of A at X, J(X;A): X -» X, a single valued, 
everywhere defined, nonlinear operator on X can be defined as J(X;A) = 
(I + XA)' 1 . 
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A two parameter family of nonlinear operators (U(t,s): 0 < s S t < T) 
defined on a subset Cl C X is called a nonlinear evolution system on Cl if for 
each x € Cl we have U(t,s)x G Cl, U(s,s)x = x and U(t,r)U(r,s)x = U(t,s)x for 
OSsSrSt^T and U(t,s)x is continuous from the triangle A = {[s,t]: 
0 t s tt $T) into X. 

A strongly continuous function x: [0,t] -» X is called a strong solution to 
the quasiautonomous initial value problem 

(2.3) x(t) + Ax(t) 9 f(t), 0 < t $ T 

(2.4) x(0) = x° 

where f: [0,T] -* X and x° € X if x is absolutely continuous on compact 
subintervals of (0,T), differentiable almost everywhere and satisfies f(t) - 
x(t) G Ax(t) for almost every t G [0,T] and x(0) = x°. 

Theorem 2.1. Let X be a Banach space with norm |-| x and suppose that A: X - 
2 X and f: [0,T] -* X appearing in (2.3) satisfy 

(1) there exists an u> G R for which the operator A + wl is m -accretive, 

(2) f G L^TjX). 

Then a unique , nonlinear evolution system (U(t,s): 0 S s $ t S T) on Dom(A) can 
be constructed which satisfies 

(i) |U(t,s)<#> - U(t,sM x « e^ (i) * * 4 * - 8 ^ - 4H X , for <p,<l> G Dom(A) and 

0 C s (t t T, 

(ii) |U(s+t,s)</> - U(r+t,r)^| x « 2 [ e^'^I^T+s) - f(T+r)| x dT, 

•’o 

for all <p G Dom(A) and all t > 0 such that s+t, r+t?T. 
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(iii) if x° € Dom( A) and the initial value problem (2.3), (2.4) has a strong 
solution x, then 

x(t) = U(t,s)x(s), for 0 ^ s U « T. 

When x° € Dom( A) , the strongly continuous function x: [0,T] -* X given by x(t) 
U(t,0)x° is referred to as a mild or generalized solution to (2.3), (2.4). 

Theorem 2.1 is a direct consequence of results given by Crandall and 
Evans and Evans in [7] and [9]. Henceforth, we shall assume that A 0 : X 0 -» 2 x o 
and f 0 : [0,T] -*■ X, satisfy (1) and (2) in the statement of Theorem 2.1 and 

that x 0 € Dom(A 0 ). We then let (U 0 (t,s): 0 $ s S t $ T} denote the 
corresponding nonlinear evolution system on Dom(A 0 ) and consider the 
approximation of mild solutions to the initial value problem (2.1), (2.2). 

Our approximation result is in the spirit of those given for nonlinear 
semigroups and evolution systems by Crandall and Pazy in [8] and Goldstein 
in [10]. However, our theorem differs from these earlier treatments in two 
ways. First, we require that the time dependent perturbation f Q be only as 
opposed to it being continuous as in [8] and it satisfying a Lipschitz-like 
condition in [10]. This distinction is especially relevant in the case of control 
systems where discontinuous input is common. The second difference is that 
we give our result in a form that is most appropriate for application to the 
development of a general approximation theory or framework and 
computational schemes for the parameter identification problems to be 
discussed in the next section. 

We shall require some set theoretic notation. For sets H n , n=0,l,2,.„, by 
lim H n DH fl we shall mean: Given x Q eH 0 , there exist x n €H n such that x n -»x Q as 
n ®. 
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Theorcm 2.2. For each n € Z + = {1,2,3, ...} let X n be a closed linear subspace of X Q . 

x 

For n = 0,1,..., let A n : X n -* 2 n be a possibly multivalued nonlinear operator on X n , and 
let f n : [0,T ]-* X n bean X n -valued measurable function defined on [0,T]. Suppose 
that there exists an € R, independent of n, for which the operators A n + w Q I 
are m -accretive, that there existsa function g€L 1 (0,T;X 0 ) for which |f n (t)| S g(t), 
a.e.t€ [0 ,T], and that lim D n D D 0 where D n = Dom(A n ) and D Q = Dom(A 0 ). 
Suppose further that for some X Q > 0 we have 

(2.5) lim J(X 0 ; A„ + w 0 I )*„ = J(X Q ; A 0 + u Q I )4> 0 

n 

whenever <p € X with lim <p„ = &.€ X. , 

n n n U U 

n 

and that 

lim f n (t) = f 0 (t) for a.e.te[0,T]. 

n -K® 

Then for each n € Z + there exists a unique nonlinear evolution system (U n (t,s): 
0 $ s ^ t $ T} on D n corresponding (in the sense of Theorem 2.1) to A n and f n 
and for 6 € with lim = 4> n G D n we have 

n n n -+oo n u u 

(2.6) lim U n (t,s)<#> n = U o (t,s)0 o , OCsUa, 

n"*°° 

with the limit being uniform in t for t G [s,T]. 

Proof. We follow Goldstein (see [10], [11]) and use an approach first suggested 
by Kisynski [13] for demonstrating the convergence of approximations to 
linear semigroups, to prove the theorem via an application of our existence 


result. Theorem 2.1. 
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Let x = {x = {x n }” =0 : x n e X n , n = 0,1,2 and Hm x n = x Q } and for 

x € x set llxll = sup|xj 0 . Then II • II defines a norm on the linear vector space 

n 

X, and the space X together with the norm II - U is a Banach space. Define the 
operator A: X - 2 X by 

dom(/l) = {x = {xX =0 ex: x n € Dom(A n ), and for each n = 1,2, ... 
there exists a y n € A n x n such that lim y n = y 0 € A 0 x 0 ), 

for X e Dom(/l), y = {y n }~=0 € A * if and 0nly if y n e A n X n’ 
n = 0,1,2, ... and lim y n = y Q . 

n - * 00 


Define an essentially x-valued function f on the interval [0,T] by f(t) = 
{f n (t)}“_ 0 . The assumptions on the f n are such that f n (t) -* f Q (t) for almost 
every t e [0,T]. However, by appropriately redefining on a set of measure 


zero, we may infer from the assumptions on the functions f that f:[0,T] -* x 
with f € L x (0,T;x). 

It is readily seen that the operator A + w 0 I is m-accretive. Let x 1 = 

x2 = < x2 Co e Domoo and let y 1 = ( y X=o e A *' and y 2 = {y 2 >r=o e 

Ax 2 . Since for each n = 0,1,2, ..., A n + u) Q I is assumed to be m-accretive, for 
X > 0 we have 


Wx 1 - x 2 P = sup |x* - x 2 | 0 S sup|x* - x 2 + X(yJ; + u 0 x* - (y 2 + U 0 x 2 ))| 0 


A a A i A | 

x + \(y + U) Q x 


Jo x 2 ))ll, 


and therefore that A + u> 0 I is accretive. Now let y = {y n }”_ 0 € x and set x = 
( x nC=o with X n = J ( x 0’ A n + n = °’ 1 ’ 2 ’ - where * 0 is chosen as in (2.5). It 

is immediately clear that for each n = 0,1,2, ..., x n € Dom(A n ) C X n . Since 
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y € x we have lim y = y ft and therefore, by assumption (2.5), that lim x„ = x n 

n _+oo n u n“*°° n U 

or x € x . Setting z n = (y n - (1 + ^ 0 w 0 )x n )/X 0 , n = 0,1,2, ..., it follows that 

A 

z e A x and lim z = z. € A n x n . We conclude that x e Dom(/l), (I + 

n n n n u u u 

n 

X 0 (^+w 0 I))x B y, and that R(I + \ Q (A + 0 ) Q I)) = x 

We have shown that the operator A and the function f satisfy conditions 
(1) and (2) given in the statement of Theorem 2.1. Therefore, a unique 

nonlinear evolution system (f/(t,s): 0 S s * t * T} on Dom(/4) corresponding to 
A and f can be constructed with U( t,s) = (U n ( t,s)}~_ 0 . Using assumption (2.5) 
it can be shown that Dom(^) = {x = {* n }“_ 0 € x : x n € D n , n = 0,1,2, ... and 
lim x = x ft ). Since R(C/(t,s)) C x , it follows that 

n-*® 

(2.7) lim U n (t,s)0 n = u 0 ( t,5)« 0 , OtstMT 

n* 4 ® 

whenever <t> € D and lim 4> = <t> n € Since each of the operators A and 

n n —+oo n 0 0 n 

n 

the functions f n satisfy conditions (1) and (2) of Theorem 2.1, unique 

nonlinear evolution systems (U n (t,s): M s It ( T) on D n corresponding to A n 

00 

and f n can be constructed. Recalling that Dom(-4) c x D n , we may 
define the family of operators (V(t,s): 0 ^ s U ( T) on Dom(/i) by 

(2.8) V(t,s)x = {V n (t,s)x n }- =0 = {U n (t,s)x„}” =0 

for x = {x n }“ =0 € Dom(yi). Uniqueness (see [9]) dictates that for each n = 
0,1,2, ..., U n (t,s)x n = V n (t,s)x n whenever (x n }~ =0 € Dom(/l). This together 
with (2.7) and (2.8) establish (2.6). The fact that the convergence in (2.6) 
is uniform in t for t € [s,T] is argued exactly as it was for the 
convergence of approximations to nonlinear semigroups in the proof of 
Theorem 3.2 in [10], 
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We note that (2.5) is also a necessary condition for the conclusion to 


hold (see, for example. Theorem 1 in [14]). 
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3. An Approximation Theory for Identification Problems 

Let H be a real Hilbert space with inner product <•,•> and 
corresponding norm |-|. Let V be a reflexive real Banach space with 
norm II • II and let V* be its dual. (All of our theory can be developed in 
complex spaces if necessary; see [6].) We denote the usual dual norm on V* 
by II -II* and assume that V is densely and continuously embedded in H 
with |v| S MlIvU, v € V, for some positive constant u. Identifying H with 

its dual, we obtain V C H = H* C V*. For <f> € V* and v € V the duality 

pairing between 4> and v is denoted by <0,v>. When <t> € H, its pairing 

with v e V agrees with the inner product of <P with v. It follows for u € 

H and v e V that Bull* * jt|u| and ivB* « /x 2 llvB. Let Q and Z be metric 
spaces and let Q be a nonempty, sequentially compact, subset of Q. The 
spaces Q and Z, and the set Q are referred to as the parameter space, the 
observation space, and the admissible parameter set respectively. 

We recall that a single valued operator A:V-.V* is hemicontinuous if 
lim A(u+tv) = Au for all u,v€V where the limit is taken in the weak sense. 

t— +0 

For each q e Q let A(q): V -» V* be a single valued, hemicontinuous, (in 
general, nonlinear) operator satisfying: 

(A) (Continuity): For each v € V, the map q -♦ A(q)v is continuous from 
Q C Q into V*. 

(B) (Equi V- monotonicity): There exist an u € R and an a > 0, both 
independent of q e Q, such that 

<A(q)u - A(q)v, u-v> + u|u-w| 2 £ aBu— vll 2 . 


for every u,v e V. 
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(C) (Equi-boundedness): There exist a constant 0 > 0, independent of 

q € Q such that 

II A(q)v U. S 0( II v H + 1), 
for every v € V. 

For each q e Q, let f( ;q) € L 1 (0,T;H) and u°(q) e H and assume that 

the mapping q -* u°(q) is continuous from Q C Q into H and that the 

mapping q -* f(t;q) is continuous from Q C Q into H for almost every t € 
[0,T], Also, for every z e Z, let u -* <t>(u;z) be a continuous map from 

C(0,T;H) into R + . 

We consider parameter identification or inverse problems of the form: 

(ID) Given observations z € Z, determine parameters q € Q which 

minimize 

<Kq) = <Ku 0 (q);z) 

where u Q (q) = u Q ( ;q) is a mild solution to the initial value problem 

(3.1) u(t) + A(q)u(t) = f(t;q), 0 < t t T, 

(3.2) u(0) = u°(q) 

corresponding to q € Q. 

By a mild solution to (3.1), (3.2) we mean a solution in the sense of 
Theorem 2.1. To be more precise, for each q € Q we define the operator 
A Q (q): Dom(A 0 (q)) C H -* H to be the restriction of the operator A(q) to the 
subset of V given by Dom(A 0 (q)) = {v € V: A(q)v e H), and prove the 

following theorem. 
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Thcorem 3.1. For each q € Q the operator A 0 (q): Dom(A Q (q)) C H -* H is 
densely defined and the operator A Q (q) + ul is m-accretive. 

Proof. We first show that for each q € Q the operator A(q) + wl: V -* V* is 
coercive. If (v } C V with lim it v U = °®, then from assumptions (B) and (C) we 

n n~ »oo n 

obtain 

lim <(A(q) + wI)v n ,v n >/Hv n ll 

n "* 00 

= lim {(<A(q)v n - A(q)0,v n > + ^vj 2 )/ H v„fi + <A(q)0,v n >/llv n B} 
n -*00 

> lim {otllv n ll 2 / Ilv n l1 - !<A(q)0,v n >|/ B v n B} 

£ lim (allv II - II A(q)9 II *} > lim aBv 1 — 0 = ® 

400 ^ n - 4 C 0 ^ 

where 0 denotes the zero vector in V. It follows that for each X > 0, the 
operator I + X(A(q) + ul): V -* V* is monotone, everywhere defined on V, 
hemicontinuous, and coercive. Consequently IR(I + X(A(q) + wl)) = V* (see 
Barbu [6], Theorem II. 1.3) and therefore R(I + X(A 0 (q) + wl)) = H. Also, for 
u,v e Dom(A 0 (q)), we may use assumption (B) to conclude 

^1 + — j|u— v| 2 $ |u-v| 2 + X oc II u — v II 2 

S |u-v| 2 + \<(A(q) + ul)u - (A(q) + ul)v, u-v> 

= <(I + X(A(q) + ul))u - (I + X(A(q) + wl))v,u-v> 

« 1(1 + X(A 0 (q) + wl)u - (I + X(A 0 (q) + ul)v| |u-v| 
or 

I u v | S |u - V + X((A 0 (q) + ul)u - (A Q (q) + ul)v)| 
which proves that A Q (q) + wl is m-accretive on Dom(A Q (q)) C H. 


To show Dom( A Q ( q) ) = H, we let u € H and for each n = 1,2, ... we set 
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u n = J(l/n;A 0 (q)+uI)u € Dom(A 0 (q)). Then, arguing as we have above, we 
find 

|u n | 2 + (l/n)a«u n i 2 $ <u - (l/n)A(q)0,u n > 

« |u| |u n | + (l/n)IA(q)ei.lu n l 

where 0 is again the zero vector in V. But then 

(3.3) (l/2)|uj* + (l/n)(a/2)lu n l 2 « (l/2)|u| 2 + (l/n)(l/2«)lA(q)0l5 

* (l/2)|u| J + (l/n)(3 2 /2ot), 

from which it immediately follows that the u n are uniformly bounded in H. 
Indeed, from (3.3) we see that (l/n)Iu n l 2 and, hence lu n l/iTn, is bounded so 
that iu l/n-*0 as n-°°. 

Also, assumption (C) yields 

Hu n - ul. = (l/n)#(A 0 (q) + ul)u n l, « (l/n){(3 + WM 2 )Hu n ! + 6). 

Since the last term in the estimate above tends to zero as n -* ®, we find 
u n -* u in V* as n -» This, together with the fact that V is dense in H 
imply that u -♦ u weakly in H as n -» 40 from which Dom(A Q (q))=H 
immediately follows. 

In light of Theorem 3.1, we may apply Theorem 2.1 with X = H, A = 
A 0 (q) and f = f( ;q). We conclude that there exist a unique nonlinear 
evolution system (U 0 (t,s;q): OSsStSTJonH satisfying (i), (ii) and (iii). 
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The mild solution u 0 ( ;q): [0,T] -»H to the initial value problem (3.1), (3.2) is given 
by u 0 (t;q) = U o (t,0;q)u°(q) for t e[0,T]. 

Remark. Under additional hypotheses on f(-;q) and u°(q) other existence 
results can be applied to obtain somewhat different notions of a solution to 
the initial value problem (3.1), (3.2). For example (see [6, p.140-144]) if 
f( ;q) € W 1,1 (0,T;H) and u°(q) € Dom(A Q (q)), then there exists a unique u(-,q): 
[0,T] - V satisfying u(-;q) € W^O.TjH), A(q)u(-;q) € LJ0,T;H) and u(t;q) + 
A(q)u(t;q) = f(t;q) a.e. t € [0,TJ. Or, if u°(q) € H and f ( - ;q) € L 2 (0,T;V*) 
then there exists a unique u( ;q) which is V*-valued absolutely continuous 
almost everywhere on [0,T], u( ;q) € C(0,T;H) n L 2 (0,T;V), u(-;q) € L 2 (0,T;V*) 
and u(t;q) + A(q)u(t;q) = f(t;q), a.e. t e [0,T], If, in addition, the mapping t -* 
t^f’(t;q) is an element in L 2 (0,T;V*) for some y > 1, then the mapping t -* 
Pu(t;q) is in L 2 (0,T;V) n L^O.T;!!). In particular, when f( ;q) = 0, the 
nonlinear semigroup (S 0 (t;q): 0 tt tT) on H defined by S Q (t;q) = U Q (t;0;q), 
t € [0,T], with generator -A Q (q) behaves like a holomorphic linear semigroup 
in that it smooths. That is, S 0 (t;q)u°(q) € Dom(A 0 (q)), t € (0,T], and the 
mapping t -* t gj- S(t;q)u°(q) is an element in L^O.TjH) for every u°(q) € H. 
Also, some generalizations are possible. For example, in assumption (B). the 
term allu - v II 2 can be replaced by a term of the form oc( II u — v II ) II u - vll where 
cc() is a continuous, strictly increasing function on 
[0,®) satisfying a(0) = 0 and lim o(x) = Or, the terms Uu - vll 2 in (B) 

x -* 00 

and II vll in (C) can be replaced by llu - v II p and lvll p '\ respectively, for 
any p £ 2. 
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The development of computational methods for the solution of the 
i infinite dimensional optimization problem (ID) requires the finite dimensional 

approximation of the abstract initial value problem (3.1), (3.2). The general 
framework that we are proposing is based upon a classical Galerkin approach. 
’ For each n = 1,2, ... let H n denote a finite dimensional subspace of H which is 

a subset of V. Let P n : H -* H n denote the orthogonal projection of H onto H n 
with respect to the <-,•> inner product. We assume that the approximating 
subspaces H n , and the projections P n satisfy 

(D) For each v € V, lim BP v — vll = 0. 

n-*“ n 

Note that assumption (D) and V densely and continuously embedded in H 
imply that lim |P u — u| = 0 for each u € H. 

n-*~ n 

For each q € Q and n = 1,2, ... we define the single-valued operator A n (q): 
H -* H by A (q)u = v for u € H where v satisfies 

n n J n n n n n 

<A(q)u ,w > = <v ,w >, w € H . 

That A n (q) is a well defined operator from H n into H n follows from the Riesz 
Representation Theorem applied to the Hilbert space H„ and the bounded 
linear functional <A(q)u n , •> on H n . Also, define f n (;q): [0,T] - H n and 
u°(q) € H n by f n (t;q) = P„f(t;q), 0 « t $ T, and u°(q) = P n u°(q), respectively. 
Note that f n (-;q) € L x (0,T;H n ) C L x (0,T;H) and that |f„(t;q)| «|f(t;q)| for q €Q 
and almost every te [0,T], 

We consider the sequence of approximating identification problems given 


by: 
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(ID n ) Given observations z € Z, determine parameters q n € Q which 
minimize 

<*> n (q) = <Ku n (q);z) 

where u n (q) = u n ( ;q) is a mild solution to the initial value problem 
in H 

n 

(3.4) u n (t) + A n (q)u n (t) = f n (t;q), 0 < t S T 

(3.5) u n (0) = u°(q) 
corresponding to q € Q. 

From the definition of the A n (q) and the assumptions (B) and (C) on A(q), 
using arguments analogous to those used to prove Theorem 3.1, it can be 
shown that the operators A n (q) + ul are m-accretive on H n . It then follows 
from Theorem 2.1 that for each n = 1,2, ... there exists a unique nonlinear 
evolution system (U n (t,s;q): 0 1 s (t IT) on H n satisfying (i) - (iii) in the 
statement of that theorem with X = H n , f(t) = f n (t;q), and x° = u°(q). The 
mild solution to the initial value problem (3.4), (3.5) is given by u n (t;q) = 
Un(t,°;q)u° n (q), t € [0,T], 

If we assume for the moment that the approximating identification 
problems (ID n ) have solutions q n € Q, then it is desirable that they in some 
sense approximate a solution q to the original identification problem (ID). 
This is in fact the case. For suppose that it can be shown that for any 
sequence (q n ) C Q with lim q n = q € Q we have 

(3.6) limu n (q n ) = u 0 (q 0 ) in C(0,T;H). 
n -*<» 


Then {q n } C Q and Q a compact subset of the metric space Q imply that there 
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exist a subsequence (q n } C {qj and a q € Q such that lim q » q. For any 

j n j-* 00 n j 

q G Q the continuity of $ implies 


<Kq) = <Ku 0 (q);z) = d> ( lim u (q );z) 

j-*® j j 


= j*™ *( u „.( q „« );z) = K.K) 


* lim <P (q) = lim <J<u (q);z) 

j — >CO 11 j j -400 11 j 

= * (lim u (q);z) = <J<u 0 (q);z) 

j -400 4, j v 

= <Kq). 


Note that in the discussion above we did not assume that a solution to 
problem (ID) exists. But rather we have shown that the existence of solutions 
q n to the approximating problems (ID n ) and (3.6) imply the existence of a 
solution q to problem (ID). When the solution to problem (ID) is unique, the 
sequence {q n } itself converges to q. 

The existence of a solution q n to problem (ID n ) for each n = 1,2, ... will 
follow from the compactness of Q and the continuity of 4> once the continuous 
dependence result: lim u n (q m ) = u n (q 0 ) in C(0,T;H n ) whenever (q m ) C Q with 
= <* 0 > has been established. Although continuous dependence for the 

m - *°° 

finite dimensional systems (3.4), (3.5) could be demonstrated via a 
modification to any one of a number of familiar continuous dependence 
results for ordinary differential equations (see, for example. Hale [12], 
Theorem 1.3.4), it is also easily handled with the approximation theory 
developed in the previous section. This and the convergence in (3.6) are 
addressed in the following theorem. 
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Theorem 3.2. If assumptions (A) - (D) hold , then 

(a) If {q n } C Q with Hm q n = q Q then lim u n (q n ) = u Q (q 0 ) in C(0,T;H), and 

(b) If {q m } C Q with lim q m = q Q then lim u n (q m ) = u n (q Q ) for each n € Z + . 

Proof. Assumption (D) and the continuity of the map q -* u°(q) from Q into H 
imply lim u°(q ) = u°(q 0 ) in H. Hence, we will have verified (a) if we can 
show that lim U n (t,s;q n )w n = U 0 (t,s;q 0 )w 0 , 0 1 s It 1 T, uniformly in t for 

n"* 00 

t G [s,T] whenever w n G K n with lim w n = w o € H. We argue this using 

n~*°° 

Theorem 2.2. Note that assumption (D) implies lim H D H and assumption 

n-*® n 

(D) together with the assumed continuity of the map q -* f(t;q) from QcQ 
into H for almost every t e[0 ,T] imply lim f„(t;q n ) = f(t;q 0 ) in H for almost 

n -frOO 

every t€ [0,T] with the f n (‘;q n ) dominated by a function geL^O.TjH) 
which is independent of n. Thus, we need only to demonstrate 
that for some X Q > 0 we have 

(3.7) nmJ(X 0 ;A n (q n ) + ul)w n = J(X 0 ;A 0 (q 0 ) + wl)w 0 
in H whenever w € H , n € Z + with lim w„ = w n . 

n n’ n -*oo n 0 

Let X 0 > 0 and set v n = J(X 0 ; A n (q n ) + wl)w n and v Q = J(X 0 ;A Q (q) + wI)w Q . 
We first show that II v n II is uniformly bounded in n. From assumption (B) we 
obtain 

X 0 allV n“ 2 « + X 0 <A ( £ ln) V n " A ^n> 9 ’ V n > 

= <d + V A n (Q„) + «I))Vn.v n >-|v n | a 
+ X o < A (q o )0 - A(q n )0,v n > - X o <A(q o )0,v n > 

= <W n’ V n > “ I V f + X 0 <A (%) e “ A ( C »n) 0 ’ V n > “ X O <A K) 0 ’ V n > 

< lw n IJv n l + X o »A(q o )0 - A(q n )eijlv n ll + X Q » A(q o )0 1 .IvJ 
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where 8 denotes the zero vector in V. This estimate together with assumption 
(C) yields 


llv „ U S ( X o 0t )’ 1 ^ W n l + a " 1,A ta n ) 0 “ A (q 0 ) 9 "* + a ' le - 

Recalling assumption (A) and that lim w = w n in H, we find that the desired 

n - * 00 n U 

uniform bound on II v n II has been established. 

Once again, from assumption (B), we find 

X O allV n- V o" 2 « X 0 W l V n“ V ol 2 + X 0 <A K) V n “ A (<O V 0’ V » " V 0> 

= X 0 U t V n" V ol 2 + X 0 <A (9 n ) v n " A (%> V 0’ V n " P „ V 0 > 

+ X 0 <A (d n ) v n - A (d 0 ) v 0 * P n V 0 " V 

+ X 0 <A(q 0 )v 0 - A(q n )v 0 ,v n - v 0 > 

= X 0 W<P n V 0 “ V 0’ V „ - V + <( l + V A .K> + U)I » V n 
- (I + X 0 (A 0 (q 0 ) + u>I))v 0 ,v n - P n v 0 > + <v 0 -v n ,v n -P n v 0 > 

+ X 0 <A(q n )v n - A(q 0 )v 0 ,P n v 0 - v 0 > 

+ X 0 <A(q 0 )v 0 - A(q n )v 0 ,v n - v Q > • 

= X 0 U<P „ V 0 “ V 0’ V n - V 0> + <W n " W 0’ V n " P „ V 0> "I V n _P n V 0 I 2 

+ " A fao) V 0’ P n V 0 " V 0> 

+ X 0 <A(q 0 )v 0 - A(q n )v 0 ,v n - v 0 > 

< X 0 ulP n V 0" V 0 l * lv n“ V 0 1 + 1 W n^ V 0 '* 1 

+ lw n -w 0 l,IP B v 0 -v 0 l + X 0 HA(q n )v n - A(q 0 )v 0 MP n v 0 - v 0 « 

+ X 0 IIA(q 0 )v 0 - A(q n )v 0 Mv n - v 0 H. 


l , n , 

The estimate ab 5 — a + — b for any h > 0 and assumption (C) allow us to 
2n 2 


argue 


X Of 3^2 y O 

— 0 — II v — v II 2 $ °-lP v -v.llj + llw -w + II w — w n H * II P v —v II 

2 n 0 n 0 0 * 2\cc n 0 n 0 nOO 
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+ X 0 IIA(q n )v n - A(q 0 )v 0 MP n v 0 - v Q B + ^ Q -«A(q 0 )v 0 - A(q n )v Q B 2 


f3u 2 X 0 n 4 <X +1) f 3n 2 n 2 3 

l 2a 2 J n 0 0 L 2X a 2 ]' n °' 


2 „2 
0 C 


+ + “ V 0») + 2 > 1P n V 0 “ V 0° + ^ nA (%) V 0 “ A^n^O 1 * ‘ 


From this, the uniform bound on llv 0, lim w = w. in H and assumptions 

n „-.<*> n O 

(A) and (D) allow us to conclude lim v„ = v. in V and that (3.7) holds. 

n-*« n 0 

An analogous, but somewhat simpler argument can be used to verify (b). 
We use Theorem 2.2 to show that for n € Z + fixed, lim U (t,s:q )w = 

>oo n ’ m m 

m 

U n (t,s;q 0 )w 0 , 0 ( s U t T, uniformly in t for t € [s,T] whenever w m ,w 0 € H n 
with lim i w = w 0 in H. Clearly lim f n (-;q m ) = f„(-;q 0 ) in L.(0,T;H n ) so that 

m -400 m -4 ® 

we need only to show that for some X 0 > 0, 


+ Wl > W m = J ^0’A„(q 0 ) + Ul)w 0 

m 

in H whenever lim w = w in H. Let = J(X n ;A„(q m ) + ul)w_ and v n = 

—too mo m v O’ 7 m 0 

m 

J(X 0 ;A n (q 0 ) + ul)w 0 . Then from assumption (B) 

X O a " V m _V o“ 2 5 X 0^ V m~^ol 2 + X 0 <A ^m) V m " A ^ m ) V 0> V m “ V 0 > 

= <(I + X 0 (A n (q m ) + ul))v m - (I + X 0 (A n (q 0 ) + ul))v 0 ,v m - v Q > 

~ ^nT'o^ + X 0 <A ( £ 1o) V 0 “ A ( £ l m ) V 0’ V m " V 

= <W m~ W 0’ V m _ V 0 > " l V m _V o| 2 + X 0 <A (%) V 0 “ A <<O V 0’ V m “ V 
< «w m -w 0 ljlv m -v 0 ll + X 0 BA(q 0 )v 0 - A(q m )v 0 l.lv m - v 0 fl 

or 

» V m " V o" * 7T |W m " W ol + 7 lA («lo ) V 0 “ 

AqOc a 


I 

! 
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Assumption (A) and lim w = w. in H yield the desired result and the 
theorem is proved. 

Remark. In practice, the approximating identification problems (ID n ) are solved 
using standard iterative search techniques (for example, steepest descent, Newton’s 
method, etc.) requiring the evaluation of </> n (q) for q e Q at each step. This in turn 
requires the integration of the finite dimensional initial value problem (3.4), (3.5). 
Once a basis for H n has been chosen, the solution to (3.4), (3.5) can be computed using 
any standard numerical integrator for ordinary differential systems. Also, the 
parameter space Q and the admissible parameter set Q are frequently functional in 
nature and infinite dimensional. When this is the case the set Q must also be 
discretized. Suppose that for each m = 1,2,..., I" 1 : Q CQ -Q is a continuous map with 
finite dimensional range and that lim I m (q) = q with the convergence uniform 

m-* 0 ® 

in q for q € Q. Set Q™ = I m (Q) (note that Q™ is a compact subset of Q) and 
consider the identification problems (ID™) defined to be the problems (ID n ) 

with Q replaced by Q m . It is clear that each of these problems admit a 

solution q^ and it is not difficult to argue that there exists a subsequence 

(q I "- i ) C {q™} with lim q"^® q, q a solution to problem (ID) (see, for example, 
n k " j,k-"» n k 

[4]). Once bases for H n and the range of I m have been chosen, problem (ID™) 
involves the minimization of a functional over a compact subset of Euclidean 
space subject to finite dimensional constraints. 

Remark (Nonautonomous systems). Theorems 2.1 and 2.2 remain valid for 
certain classes of temporally inhomogeneous or time dependent operators A = 
A(t). To be more precise, the family of operators A(t): X -* 2 X must be 
m-accretive on X for almost every t € [0,T] and must satisfy 
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(3.8) |J(X;A(t))x - J(X;A(s))x| x S X|h(t) - h(s)| x L(|x| x ) 

for each x € X, every X satisfying 0 < X $ X Q for some X Q > 0, some h € 
L 1 (0,T;X), some continuous, non decreasing function L: [0,®) - [0,°°) and 
almost every t,s € [0,T] (see [8], [9]). (Note that for simplicity we have taken 
w = 0; however, the discussion to follow remains valid for any u € R.) The 
primary motivation for developing the framework outlined above was to 
define readily verifiable conditions on the operators A(q): V -* V* that if 
satisfied would (i) also automatically be satisfied by the Galerkin 
approximation A n (q) and (ii) lead to the desired convergence of solutions to 
the approximating identification problems to a solution to problem (ID). The 
natural assumption to add to (A) - (C) that certainly satisfies criterion (i) and 
that could conceivably lead to an estimate of the form (3.8) in H is that 

(3.9) IA(t;q)v - A(s;q)vH, $ |h(t) - h(s)|L(|v|) 

for each v € V, almost every s,t € [0,T] and some h € Lj(0,T;H) and some 
continuous nondecreasing L: [0,®) -* [0,®), both of which do not depend upon 
q € Q. Unfortunately, however, we can only show that (3.9) leads to an 
estimate of the form 

(3.10) |J(X;A 0 (t;q))u - J(X;A 0 (s;q))u| « ✓! |h(t) - h(s)| L(|u|) 

for each u € H. Moreover, it is not clear to us how, or if, the proof of the 
fundamental Theorem 2.1 given in [9] could be modified so that (3.10) would 
suffice. We have explored alternative approaches and developed other 
techniques for treating the nonautonomous case (for example, in the linear 
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case, based upon some ideas in Tanabe [18], and in the strongly monotone case, 
via a variational formulation which can be found in Barbu [6]). These 
results will appear soon in forthcoming papers. 
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4. Applications and Examples 

We briefly describe some classes of systems to which the general 
framework developed in the previous section applies. In our discussion below 
we consider theoretical aspects only. Implementation questions will be treated 
and the results of our numerical studies will be reported on elsewhere. 

Example 4.1. Linear regularly dissipative operators. The approximation 
theory for inverse problems for systems involving linear regularly dissipative 
operators was treated in detail by Banks and Ito in, and is the central focus 
of, [2] and [3]. We show here that the linear theory is a special case of the 
nonlinear theory given in Section 3. 

Let the spaces H, V, V* and Q and the set Q be as they have been defined 
above. For each q € Q let a(q)(-,*) be a sesquilinear form defined on V x V 
which satisfies the conditions: 

(A*) For each v e V the mapping q -* a(q)(*,v) is continuous from Q C Q 
into V*. That is given € > 0 there exists a 6 > 0 such that 

sup |a(q 0 )(u,v) - a(q)(u,v)| < e 

u€V u 

Hu 11=1 

whenever d(q Q ,q) < 5 where d denotes the metric on Q. 

(B*) There exist an u € R and an a > 0, both independent of q € Q, for 
which a(q)(v,v) + w|v| 2 £ aBvB 2 for every v e V. 

(C 1 ) There exists a constant B > 0, independent of q G Q, such that 

|a(q)(u,v)| S BlluH Ivl for every u,v G V. 


-25- 


When conditions (A 1 ) - (C 1 ) are satisfied it is not difficult to argue that for 
each q € Q an operator A(q) € £(V,V*) can be defined by 

[A(q)v](u) = <A(q)v,u > = a(q)(u,v), 

u,v e V and that A(q): V -» V* satisfies (A) - (C). It then follows from 
Theorems 2.1 and 3.1 that there exists a unique nonlinear evolution system 
{U Q (t,s;q): 0 ^ s U ^ T) on H corresponding to the initial value problem 

u(t) + A 0 (q)u(t) = f(t;q), 0 < t S T 
u(0) = u°(q) 

where for each q € Q, f(-;q) € L^O.TjH), u°(q) € H and A 0 (q): Dom(A 0 (q)) C 
H -* H is the restriction of A(q) to the set Dom(A 0 (q)) = {v € V: A(q)v € H}. 
The operator — A 0 (q) is the infinitesimal generator of an analytic semigroup 
(T 0 (t;q): t > 0} on H (see [18]) and for <t> € H 

(4.1) U 0 (t,s;q)4> = T 0 (t-s;q)tf> + f T 0 (t-T;q)f(T;q)dT. 

It can be shown that the semigroup (T 0 (t;q): t > 0} admits an extension (T(t;q): 
t £ 0} which is an analytic semigroup on V* with generator A(q): V C V* -» V*. 
Also the restriction of (T 0 (t;q): t £ 0} to V, call it (T(t;q): t > 0), is an analytic 

»v <v 

semigroup on V with generator A(q): Dom(A(q)) C V -• V, the restriction of 
A(q) to the set Dom(A(q)) = (v e V: A(q)v e V) (see [3], [18]). Consequently, 
with appropriate assumptions on f( ;q), the evolution system (U Q (t,s;q): 0 « s 
« t S T) admits an extension (U(t,s;q): 0 ^ s It 1 T) which is an evolution 
system on V* and a restriction (U(t,s;q): 0 1 s It 1 T) which is an evolution 
system on V. 
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It is clear from (4.1) that when A(q) is linear, we may take f(;q) = 0 and 
consider only the approximation of the semigroup (T 0 (t;q): t £ 0). For each 
n = 1,2, ... let the finite dimensional subspaces H n of H and the corresponding 
orthogonal projections P n be as they were defined in Section 3 and assume 
that condition (D) is satisfied. Denote the Galerkin approximations to A(q) 
(i.e. the restriction of A(q) to an operator from H n into H* = H n ) by A n (q) and 
set T n (t;q) = exp(— tA n (q)), t £ 0. Theorem 3.2 then implies that 

(4.2) Urn |T n (t;q n )P n u°(q n ) - T 0 (t;q 0 )u°(q 0 )| = 0 

uniformly in t, for t € [0,T] whenever {q n } C Q with lim q n = q Q € Q, and 
the mapping q -* u c (q) is continuous from Q C Q into H. In addition, recalling 
that we required that H n C V for all n = 1,2, .... an inspection of the proof of 
Theorem 3.2 reveals that in the linear case with the existence of the semigroup 
(T(t;q): t > 0) on V, we may apply Theorem 2.2 with X = V and conclude that 

(4.3) lim lT n (t;q n )P n u°(q n ) - T (t;q 0 )u°(q 0 ) I = 0 

uniformly in t for t € [0,T] whenever lim q n = q 0 , u°(q) € V and the map 
q -* u°(q) is continuous from Q into V (see also [3]). Then for <t> € H, setting 

U n (t,s;q)P n * = T n (t-s;q)P n tf> + f T n (t-T;q)P n f(T;q)dT 

J 8 

under appropriate assumptions on f( ;q), (4.2) and (4.3) continue to hold with 
T n (t;q), T 0 (t;q), and T(t;q) replaced by U n (t,s;q), U 0 (t,s;q), and U(t,s;q), 
respectively, with the convergence being uniform in t, for t € [s,T]. Hence the 
linear theory and results of [3] are a special case of the nonlinear theory of 


Section 3. 
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We note that in the context of the identification problem, the fact that 

the stronger V-convcrgence given in (4.3) can be obtained is significant. 

Indeed, (4.3) permits the relaxation of the continuity assumption on the 

performance index 4> to the requirement that for each z € Z, the mapping u -» 

4>(u,z) be continuous from C([0,T];V) into R + . This can have the effect of 

significantly enlarging the class of allowable observations. For example, in 

the case of a one dimensional parabolic system formulated in H = L 2 with V 

in H 1 , spatially discrete (i.e. pointwise, as opposed to distributed in space) 

measurements will suffice (see [3] and (5]). 

Among the class of linear regularly dissipative operators which arise from 

a form satisfying (A 1 ) - (C') are the familiar elliptic partial differential 

# i 2 +i+i 

operators on L,. Briefly, let ft be a region in R and let Q - * L w (ft). Let 

1 n=l 

Q be a compact subset of Q with the property that if q= {( a U’ b ,c): iJ-l_l)€Q, 
then for some a > 0 independent of q € Q. 
i 

£ a ij(x)C i C; > «UI 2 

i,j=l 1J 1 J 


for every x € ft, and every £ € R*. For q € Q and u,v € H x (ft) set 


a(q)(u,v) 


»(*) 


b ^> ax 


9u(x) 9v(x) 


s 

+ I 


9xj 
9u(x) 


9x. 


i=l 


v(x) + c(x)u(x)v(x)}dx 


with H = L 2 (ft) and V any closed subspace of H x (ft) containing Hj(ft), it can be 
shown (see [18]) that a(q)(-, ) satisfies (A’) - (C 1 ). The operator A(q) is given 
formally by 
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* 9 9 1 9 

(4.4) A(q) = - Z - — a;.(x) - — + E b;(x) - — + c(x). 

i,j=l 0Xj 1J OX . i=l 1 9Xj 

When 9ft is sufficiently smooth, A(q) is the elliptic operator given by (4.4), and 

V is chosen to be either Hj(ft) or H x (ft), the equation (3.1) becomes a parabolic 
partial differential equation with either Dirichlet or Neumann boundary 
conditions. 

For H = L 2 (ft) and V a subspace of H^ft), choosing the approximating 
subspaces to be the span of an appropriate collection of first order spline 
functions will typically satisfy assumption (D) (see [15] and Example 4.2 
below). 

Example 4.2. Nonlinear Elliptic Operators. Let 0 be a bounded region in R* 
with smooth boundary T = 9ft. For a = (a .... cc g ) a multi-index, let lot) = ct, + 
a 2 + • • - + aj and denote the otth order generalized, or distributional 

derivative of a function u by D'Sj; that is, 
a, a. 

„ 9 1 9 

D a u(x) = — — • ■ • — — u(x), x € ft . 
ax, 1 9x, J 

Let m be a nonnegative integer and let 8u denote the vector valued function 
of length N = (*J m ) whose components are all of the partial derivatives of u 
of order greater than or equal to zero and less than or equal to m. 

For each multi-index a with |a| « m, let (x,0 -♦ a a (x,0 be a real valued 
function defined on ft x R N which is measurable in x and continuous in £. We 
assume that 

(1) there exist age L 2 (ft) and a positive constant 7 such that 
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(4.5) |a a (x,OI * 7(151 + g(x)) 

for almost every x e ft, each £ € R N and all a with |a| « m, and 
(2) there exists a positive constant \ such that 

(4 - 6) ( a «( x .5) - a (x,h))(5 - o ) > X E I5 a - hj 2 

|a|^m |a|^m ^ 

for almost every x e ft and all (;,n € R N . 

Let H = L 2 (ft) and let V be any closed subspace of H m (ft) which 
contains H™(ft). Define the operator A: V -» V* by 

(4.7) (Au)(v) = I [ a (x,5u(x))D ot v(x)dx, 

for u,v € V. The operator A given by (4.7) is the distributional form of the 
formal differential operator 

(4.8) (Au)(x) = E (-l) a D 0t a <x (x,6u(x)). 

|a|$m 

A differential operator of the form (4.8) is referred to as a nonlinear elliptic 
operator and the partial differential equation 

(4.9) ^ (t,x) + E (-l) a D a a (x,6u(t,x)) = f(t,x) 

at |a|$m 

is said to be of nonlinear parabolic type. When V = a solution in V* 

to the abstract equation 

u(t) + Au(t) = f(t) 

with A given by (4.7) corresponds to a variational solution to (4.9) which 
satisfies Dirichlct boundary conditions. When V = H m (ft), a variational 
solution to the Neumann problem is obtained. Note that in the linear case we 


I 
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have 


a (x,6u(x)) = E a a > e (x)D 6 u(x). 
a |B|«m 

Under the assumptions above, it is not difficult to show that A given by 
(4.7) is hemicontinuous and satisfies conditions (B) and (C) given in Section 3. 
With an appropriate choice of the space Q and the set Q, condition (A) can be 
satisfied as well. 

A quasilinear model for heat conduction or mass transfer in which the 
heat or mass flux is a function of the temperature or mass fraction gradient 
discussed in [16] and [17] leads to a nonlinear elliptic operator and a nonlinear 
parabolic partial differential equation of the forms (4.8) and (4.9), respectively 
with m = 1. Let Q be a bounded region in R* with smooth boundary and let 
Q = Lctffl x R*). Let Q be a compact subset of Q with the property that 
q € Q if and only if the mapping £ -* q(x,0 is C 1 on R* for almost every 
x € n and there exists a X > 0 (which does not depend on q) such that 

(4.10) e.^q(x,0| 5=e -a - h) + q(x,0)U. - n.) * X(C. - n,), 

for i = 1,2, .... I, almost every x € O and all 0,5, n € R*. (When 1 = 1, the 
function q(x,$) = q(0 = (1 - -5e'^ 2 ) satisfies (4.10).) 

Let H = Lj(O) and let V be any closed subspace of H x (0) which contains 
Hj(Q). Then V C H C V* and for each q € Q define A(q): V - V* by 

(4.11) (A(q)u)(v) = [ q(x,Vu(x))7u(x)’ W(x)dx 

J n 

for u,v € V. Note that for each q € Q the operator given by (4.11) is of the 
form (4.7) with 
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(4.12) a a (x, 6 u(x)) = q(x,Vu(x))D“u(x) 

for x € O and all a with |a| = 1 and a a = 0 for |o 4 = 0. The nonlinear 
parabolic partial differential equation (4.9) takes the form 

0U 

— (t,x) - V.q(x,Vu(t,x))Vu(t,x) = f(t,x), t > 0 , x € n . 

Taking II • II to be the usual norm on H 1 (ft), it follows that 

BA(q 0 )u - A(q 1 )ull. $ |q 0 - q^JIull 

for each u € V and q 0 ,q x € Q. Since Q is a compact subset of x R*), it is 
easily verified that a a given by (4.12) satisfies a growth condition of the form 
(4.5) with 7 and g independent of q € Q. An application of the. .mean value 
theorem together with assumption (4.10) imply the existence of a X > 0, 
independent of q €Q, for which (4.6) holds. Consequently the conditions (A), (B), 
and (C) given in Section 3 are satisfied, and our general theory can be applied. 

With regard to approximation, polynomial spline function based Galerkin 
subspaces can often be shown to satisfy condition (D). For example, when 
i = 1 and H = ( 0 , 1 ) in the nonlinear heat conduction/mass transfer example 
discussed above, the subspaces H n can be chosen as the span of the linear 
B-spline ("hat") functions with respect to the uniform mesh (0, 1 /n, 2/n, ..., 1 } 
appropriately modified to satisfy stable, or geometric, boundary conditions. 
Familiar error estimates for interpolation and the Schmidt inequality can then 
be used to verify that condition (D) is satisfied (see [5]). Generalization to 
higher dimensions is possible, and can often be achieved via tensor 
products of one dimensional elements (see [15]). 
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S. Concluding Remarks 

We have developed a general abstract approximation framework for the 
identification of nonlinear distributed parameter evolution systems. The class 
of systems to which our theory applies are those whose dynamics can be 
described by a nonlinear operator which satisfies conditions that are the 
natural nonlinear extensions, or analogs, of the properties of regularly 
dissipative, or abstract parabolic, linear operators. The approach we have 
taken is based upon the defining of a sequence of approximating finite 
dimensional identification problems in which the systems to be identified are 
Galerkin approximations to the original, underlying, infinite dimensional 
nonlinear dynamics. Under a weak continuity assumption with respect to the 
unknown parameters to be identified, equi-boundedness and equi- 
monotonicity conditions, and an approximation assumption on the Galerkin 
subspaces (all of which are readily verified for wide classes of nonlinear 
distributed systems and finite element subspaces), we are able to demonstrate 
that solutions to the approximating problems exist, and, in some sense, 
approximate (i.e. subsequential convergence) solutions to the original infinite 
dimensional identification problem. We have shown that the linear theory 
presented in [2] and [3] is a special case of our nonlinear framework and that 
our results are applicable to a reasonably wide class of nonlinear elliptic 
operators and corresponding nonlinear parabolic partial differential equations. 
In particular, we have considered application of our theoretical framework to 
a quasi-linear model for heat conduction or mass transport. 
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The general approximation result for nonlinear evolution systems 
discussed in Section 2 is applicable to a much broader class of nonlinear 
dynamical systems than we subsequently treated in Section 3. For example, 
this class of systems would include those with dynamics described by set 
valued maps or multifunctions, and (after minor modification to the 
general theory) time dependent or nonautonomous operators. We are 
currently investigating these features of the general approximation theory 
in the context of parameter estimation problems. Also, we would like to 
be able to weaken the somewhat restrictive strong monotonicity condition. 
Any progress that we might make in these efforts would have the potential 
to significantly enlarge the class of nonlinear systems to which our theory 
and framework would apply. Finally, extensive numerical or 

computational studies designed to demonstrate the feasibility and point out 
the limitations of our schemes and general approach are currently 
underway and will be reported on in a forthcoming paper. 
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